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Abstract We present investigations on single Ni/Si related 
color centers produced via ion implantation into single crys- 
talline type Ila CVD diamond. By testing different ion dose 
combinations we show that there is an upper limit for both 
the Ni and the Si dose (lO'^cm^^ and 10^" cm^^ resp.) due 
to creation of excess fluorescent background. We demonstrate 
creation of Ni/Si related centers showing emission in the spec- 
tral range between 767 nm and 775 nm and narrow line-widths 
of « 2nm FWHM at room temperature. Measurements of 
the intensity auto-correlation functions prove single-photon 
emission. The investigated color centers can be coarsely di- 
vided into two groups: Drawing from photon statistics and the 
degree of polarization in excitation and emission we find that 
some color centers behave as two-level, single-dipole systems 
whereas other centers exhibit three levels and contributions 
from two orthogonal dipoles. In addition, some color centers 
feature stable and bright emission with saturation count rates 
up to 78kcounts/s whereas others show fluctuating count 
rates and three-level blinking. 

PACS 42.50.Ar, 81.05.ug, 61.72.jn 



1 Introduction 

Single-photon sources are of great interest because of the 
numerous possible applications in quantum information pro- 
cessing. A prominent example is the BB84-protocol for quan- 
tum key distribution 1 1 1 that can be realized with single pho- 
tons as information carrier. For this application one requires 
a single-photon source that can be triggered, features a high 
emission rate to assure fast information transmission and ex- 
hibits a naiTow line-width to allow for an effective spectral 
signal filtering. In principle, such a source can be realized 
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by a single emitting quantum system. If such a dipole un- 
dergoes cycles of pulsed excitation, light emission and re- 
excitation before emitting the next photon, triggered single- 
photon emission can be achieved. Admittedly, many of the 
existing demonstrations e.g. single atoms (]2l or single quan- 
tum dots |3| hold experimental complexity that limits their 
suitability for practical applications. 

In comparison, color centers in diamond have larger po- 
tential for feasible single-photon sources: These point defects 
consisting of combinations of impurity atoms and lattice va- 
cancies exhibit an atom like spectrum even at room temper- 
ature (RT) |4|. Additionally the high photo-stability of many 
color centers offers the great advantage of long-time stable 
single-photon emission ISj. 

The best known color center is the so called Nitrogen- 
Vacancy center (NV) with a zero-phonon line (ZPL) at 638 
nm. Beyond single-photon emission ||6l several key experi- 
ments towards quantum information have already been per- 
formed, including e.g. the realization of a two-qubit gate Q 
and multipartite entanglement |[8| . However, NV centers have 
the detrimental property of a broad emission bandwidth of 
w lOOnm due to strong vibronic coupling to the diamond lat- 
tice. This fact strongly limits their suitability for applications 
in quantum optics and quantum information where narrow- 
band emission of single photons is essential, e.g. quantum 
key distribution or linear optics quantum computation. 

A more auspicious kind of color centers are nickel-related 
centers. NE8 centers, defects that consist of one Ni atom 
and four N atoms ||9l, show a naiTow line (1.2nm FWHM) 
around 800 nm at RT fTOl. For single NE8 centers, photon 
anti-bunching [lOJ and triggered single-photon emission LlTJI 
have been demonstrated. Investigations of color centers con- 
sisting of Ni and possibly Si in nano-crystals yielded high 
saturation count rates of 200kcounts/s at an emission wave- 
length of 768 nm with an excited-state lifetime of 2ns lfT2l . 
However, targeted placing of single emitters for applications 
in quantum information requires their deterministic creation. 
Recently the first production of single Ni/Si related color cen- 
ters via ion implantation in single crystalline diamonds has 
been reported Iil3il . We here present further investigations on 
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optimized implantation parameters of Ni/Si centers and their 
optical and structural properties. 



2 Experimental 

For our investigations we used type Ila CVD single-crystal- 
line diamonds (Element Six, Isle of Man, UK). Since the 
exact atomic composition and the creation efficiency of the 
Ni/Si related color centers are still unknown, we systemat- 
ically varied implantation energies and doses of Ni and Si 
ions to find optimal implantation parameters. We are aiming 
at a color center density that is as high as possible and still 
allows for the localization of single centers with a confocal 
microscope. We chose two different implantation energies: at 
first to investigate the influence of the implantation energy 
on the creation efficiency and secondly to explore the depen- 
dence of emission properties on the distance of color centers 
from the diamond surface: On the one hand, if the distance 
of a single emitter from the diamond surface is larger than 
approximately half the emission wavelength, refraction at the 
surface becomes relevant, lowering the collection efficiency. 
On the other hand, surface effects might act detrimental on 
color centers if the distance to the surface is too small (e.g. 
NV centers change their preferred charge state depending on 
the distance to the surface fT4lfT5]| ). An overview of the im- 
plantation parameters can be found in the following table. 



ion 


energy (keV) 


depth (nm) 


dose (ions/cm^) 


Ni 


86 


39 


1011-1014 


Si 


56 


10^-1012 


Ni 


815 


363 


1011-1014 


Si 


600 


lO^-lOi^ 



The mean implantation depths were calculated with Monte 
Carlo simulations (SRIIVU). Implantations were performed 
at the Georg- August Universitat (Gottingen, Germany) and 
at the RUBION (Ruhr-Universitat, Bochum, Germany). Af- 
ter the implantation, all samples were annealed at 1200°C in 
vacuum for one hour followed by cleaning in oxygen plasma 
to remove graphite from the sample surface formed during 
the annealing. 

Fluorescence is excited through a microscope objective 
with NA = 0.8 by a cw frequency doubled DPSS laser with 
a wavelength of 67 1 nm. A set of dielectric filters blocks the 
excitation light from entering the analyzing beam path. The 
fluorescence can be guided into a grating spectrometer or into 
a Hanbury-Brown-Twiss setup (HBT) containing two ava- 
lanche photo diodes (APDs) to measure the second order cor- 
relation function of the photoluminescence signal: 
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1 Stopping range of ions in matter: http://www.srim.org/ 



Fig. 1 Emission spectrum of a sample implanted with Ni (86keV, 
IQi-'cm^^) au(] si (56keV, IQH cm^^)^ excited with 671 nm laser 
light. Around 738 nm the emission line of the SiV center is visible, 
the feature rising till 881 nm is attributed to the so called 1.4eV cen- 
ter. Inset: detailed emission spectrum of the 1.4eV center ensemble 
measured at 30 K. 



This function represents the probability to detect two pho- 
tons separated by the time interval T. The intensity signal 
of the APDs can also be used for confocal scans of the dia- 
mond sample. Limited by the focus size of the excitation light 
a point like emitter appears with a lateral size of w 0.6 |jm 
FWHM and an axial size of 1.9|jm FWHM in the scans. The 
complete detection efficiency Tj of our setup (taking into ac- 
count the collection efficiency of the objective, the detection 
efficiency of the APDs and the transmission of all optical 
components) can be estimated as 7] =2.2%. All following 
measurements were performed at room temperature as far as 
not stated differently. 



3 Results and Discussion 

3.1 Implantation Parameters 

Photoluminescence investigations of samples implanted with 
high doses of both Ni and Si (lOi^cm^^ ^nd lOHcm^^) es- 
sentially reveal two features in the room temperature fluores- 
cence spectrum (see Fig. [T]) independent of the implantation 
energy: Around 738 nm one can observe the ZPL of an en- 
semble of SiV centers (cf. lfT6l ). Following Sternschulte et 
al. ifTTl the fluorescence up to 825 nm can be attributed to 
the SiV vibronic sideband. The feature above 825 nm with 
a maximum at w 881nm is due to so called 1.4eV centers 
iTTSJI . possibly consisting of one Ni+ ion at the center of a 
lattice di-vacancy ri9|. At cryogenic temperatures, the emis- 
sion spectrum of these color centers consists of two lines at 
883 nm and 885 nm (see inset of Fig.[T]and Ref. Il20ll ). 

In Fig.|2]the intensity of the 1 .4eV fluorescence measured 
on samples with different Ni doses is shown: The intensity 
follows a linear characteristic for the lower dose values as 
expected for a color center containing one Ni atom. It is ob- 
vious, however, that the slope reduces at the step between 
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Implantation dose (cm" ) 

Fig. 2 1 .4 eV fluorescence intensity of implanted samples in depen- 
dence on the Ni implantation dose. Inset: SiV fluorescence on a sam- 
ple homogeneously implanted with lO'^ cm^^ Si. Ni was implanted 
with four different doses in form of round circles in each comer. The 
dose decreases from lO'^cm^^ to 10^' cm^^ clockwise, starting in 
the upper left corner. 



the two highest doses 10 cm^ and 10 cm^". This behav- 
ior can be explained by increasing radiation damage influenc- 
ing the creation efficiency of color centers and/or their emis- 
sion properties negatively. The assumption of a damaged dia- 
mond lattice is supported by the observation of enhanced line 
widths of the two 1 .4eV center lines at a temperature of 30 K 
with growing implantation dose: At a dose of 10^^ cm^'^ we 
determined the widths of both lines of the doublet to be about 
0.3nm, increasing up to 0.7 nm at lO'^^cm^^. This inhomoge- 
neous broadening evidences a modified crystalline environ- 
ment for single color centers. 

The negative influence of too high Ni doses on the emis- 
sion of other color centers can be revealed by a confocal 
overview scan of one of our samples (see inset in Fig. [J]). Al- 
though the sample is homogeneously implanted with Si ions, 
a reduced SiV ZPL emission intensity in the areas of the high- 
est Ni densities is visible. These observations lead to an upper 
limit for the Ni implantation dose of 10^^ cm^^ to avoid sim- 
ilar emission quenching of single Ni/Si related centers. 

A further restriction for the ion dose is given by unwanted 
fluorescent background in the spectral range between 760 nm 
and 820 nm where we expect the Ni/Si related color center 
emission. This luminescence is produced by the vibronic fea- 
tures of the SiV centers and the fluorescence of the 1.4eV 
centers. Axial focus scans of different samples show a strong 
dependence of the unwanted fluorescence on the implanta- 
tion doses (see Fig.O - due to the strong SiV vibronic side- 
band, the effect is much stronger for SiV centers than for the 
1.4eV centers. Thus, we determine maximum implantation 
doses lO'^cm^^ for Ni and lO'^cm^^ for Si to keep the un- 
wanted fluorescence background at a tolerable level. 
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Fig. 3 Depth scans in the spectral window 765 -775 nm on eight 
different samples with different dose combinations. 



3.2 Spectroscopy & Single Photon Emission 

For the investigation of color center formation we performed 
confocal scans of our samples, detecting the fluorescence in- 
tensity in the spectral window of 765nm-820nm. Although 
we were able to observe several color centers with emission 
wavelengths varying between 786 and 815nm, the majority 
of detected color centers showed an emission wavelength be- 
tween 767 and 775 nm (examples see Fig. H)). The fact that 
these color centers can only be found in samples implanted 
with both Ni and Si (see Ref. |13| ) leads, consistent with 
Aharonovich et al. lfT2l . to the assumption that this last group 
are Ni/Si related color centers. However, the density of these 
color centers in our samples is still quite low: on average only 
one center per 2000 |jm^, measured on samples with the op 
timal doses as discussed above (lO'^cm 



-^NiandlO'^cm-^ 



Si). This corresponds to creation efficiencies of 5 x 10^^ in 
relation to the Ni dose and 5 x 10^'' in relation to the Si dose. 
Due to these low values it was not possible to analyze the 
color center density in dependence on the ion implantation 
dose. The different implantation energies and resulting depths 
of color centers did not lead to a significant difference in cre- 
ation efficiency. 

The sub-Poissonian character of the single Ni/Si centers' 
emission can be demonstrated by recording the second-order 
correlation function g^^'. In Fig. |5]the normalized intensity 
correlation measurements g^iUz) of color centers (2) and 
(4) of Fig. |4] are shown; the data have been corrected for 
background events corresponding to the signal-to-noise ra- 
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Fig. 4 Emission spectra of four Ni/Si related color centers in the 
spectral range 765 -775 nm. The spectra can be fitted with one or two 
Lorentz peak functions. The line widths of the single peaks amount 
tol.4-2.7nmFWHM. 



tios of w6:l for both emitters, a value which is obtained 

(2) 

from the confocal scans. In both cases the value of gmeas(O) 
is clearly smaller than 0.5, indicating the presence of single- 
photon emitters. For the derivation of the correct fit func- 
tions, we have to distinguish between the different emitters. 
As the correlation function of emitter (2) quickly reaches 

gJneas(f) = 1 for T > and remains flat we assume a two level 
system - in the equations below, the ground state is labeled 
as "1", the excited state as "2". Solving the rate equations for 
such a system yields 



,(2) 



(T) = l- 



M/fi 



The coefficient Ti is defined as 



Tl 



1 



'"12 + '"21 



(2) 



(3) 



with r,y representing the transition rate from state i to state j. 
Because of the timing jitter of our HBT setup we have to con- 
volute the function in Eqn.[2]with the device response func- 
tion (DRF), measured to be a Gaussian peak with a 1/Ve 

(2) 

width of 354 ps. Fitting the resulting function g^^^ to the data 
of emitter (2) yields a time constant of Ti = 1 .45 ns. 

(21 

In case of center (4), gmias exceeds the value of one for 
time scales T > 5ns. This photon bunching is a sign for the 
presence of a longer lived shelving state (signed as "3") be- 
sides ground and excited state. According to Kitson et al. 
on . the g'2) function of such a three-level system is given 
by 



jP)(t) = 1 - (1 +fl)e-l^l/^' +fle-l^l/^2 



with coefficients 



Ti,2-2/(A±VA2-4B), 

„ ^ l-'-3lT2 

'•31(^2-^1) ' 



(4) 



(5) 
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Fig. 5 Intensity correlation measurements of two color centers of 
Fig.|4]in the spectral window 765 -775 nm. The higher noise in the 
upper graph is due to the lower mean emission count rate of emitter 
(2) as compared to emitter (4). 



and the abbreviations 



A =='"12 + '"21 +'"23 + '"31, 
B = '"12'"23 + '"12'"31+'"21'"31 + '"23'"31 



(6) 



Eqn. |4] also has to be convoluted with the DRF to yield the 
correct fit function. The results for emitter (4) are Ti = 0.83 
ns, T2= 42.2ns, a = 0.16. 

(2) 

The deviations of gf^i (0) predicted by the fit functions 

(21 

and the measured values grneas(O) in both cases is smaller than 

(21 

0.027. This small value proves that the difference of gmeas(O) 
from zero is solely due to the limited timing resolution of 
our photon counting setup and the emitters show perfect anti- 
bunching. 

In summary, the measured coiTelation functions indicate 
that emitter (2) behaves as a two-level system whereas emit- 
ter (4) has to be modeled as a three-level system. A similar 
distinction of emitters has recently been observed for Cr re- 
lated color centers |22j . We note here, that emitter (2) shows 
a spectrum that can be described as a single fine, whereas the 
spectrum of emitter (4) is double-peaked with a separation 
of the two peaks of about 3.6 nm (see Fig. |4]i. The spectra 
of most investigated Si/Ni color centers fall in either one of 
these two categories. 

To gain further insight into the Si/Ni center structure we 
investigated the polarization properties of the fluorescence. 
For emitter (2) we first measured the dependence of the count 
rate on the polarization angle of the excitation light (see Fig. 
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• Emitter 2 
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30 60 90 120 150 180 

Excitation polarization angle {") 

Fig. 6 Count rates of emitters (2) and (4) in dependence on the ex- 
citation polarization. The fit functions are ordinary sines. The the 
curve of emitter (4) is shifted ~ 20° due to a tilt of the sample in our 
setup. 



|6ll. It can be well described by a sinusoidal variation; the vis- 
ibility 



y = 



(7) 



amounts to 92 %, which means, that the excitation cross sec- 
tion is nearly a perfect dipole; the deviation from 100% can 
be explained with depolarization introduced by the dichroic 
mirror in front of the objective and the objective itself. Fix- 
ing the excitation polarization to a maximum of the count 
rate and measuring the polarization degree of the light emit- 
ted from the color center yields a similar result: a periodicity 
of w 180° with maxima at sa 90° and a polarization degree of 
98.5 %. This observation is consistent with the existence of a 
single dipole. From the positions of the intensity maxima in 
both measurements and the known orientation of the fee di- 
amond lattice in our setup, one can deduce the alignment of 
the dipole or at least its projection along the (001) direction. 

On the other hand, emitter (4) shows a different behav- 
ior (see Fig. |6]i: The periodicity of the count rate in depen- 
dence on the excitation polarization likewise is about 180°. 
The deviation of the angle of minimum count rate from 90° 
can be explained with a certain tilt of the sample on the sam- 
ple holder However, the visibility amounts to 62.1 % only - a 
fact that might be attributed to the presence of a second dipole 
at an orthogonal direction. A similar behavior is known from 
the NV center [23 1 and for the 1 .4eV Ni related center Il24l . 

To summarize, one class of emitters (from now on labeled 
as "2L") behaves as two-level, single-dipole systems with a 
single emission line spectrum, whereas the other class ("3L") 
exhibits a double-peaked spectrum, a three-level system and 
the presence of two orthogonal dipoles. These observations 
are consistent with a three-level model for the Si/Ni color 
center where ground state ("1") and excited state ("2") are 
connected via a fast radiative transition, and the metastable 
shelving state ("3") lies close in energy to the excited state 
2. For emitters of class 3L the splitting of levels 2 and 3 is 
pronounced, leading to the typical three-level behavior. On 
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Fig. 7 Mean count rate of emitter (2) in time windows of 50 ms ex- 
cited with an intensity of 109kW/cm recorded for about one hour. 
Additionally to the "on" and "off" level there is an intermediate level 
at ~ 6500counts/s. 



the other hand we might assume that for the emitters 2L local 
strain in the diamond crystal or a local charge distribution 
in the vicinity of the centers shift levels 2 and 3 into close 
degeneracy, such that populations can thermalize, resulting 
in an effective two-level system. The validity of our proposed 
model has to be verified by further experiments. 

An important parameter with regard to a practical single 
photon source is the count rate of a single color center. The 
different emitters yielded different results: While the count 
rate of emitter (4) was absolutely stable, emitter (2) showed 
a fluctuating count rate. This blinking can clearly be seen 
in Fig. |7] where the count rate is shown for a time of about 
one hour in time bins of 50ms. Besides periods of maxi- 
mum count rate (~ 17000counts/s) one can see "off" peri- 
ods (Ocounts/s) and periods, where the count rate reaches an 
intermediate value of w 6500counts/s. Such a kind of three- 
level blinking has been recently observed for luminescence 
from a semiconductor quantum dot ESl . In that investigation 
the existence of two emitting and one dark level has been ex- 
plained by reversible electron trapping in surface states open- 
ing non-radiative decay channels. In case of the Ni/Si center 
an analogous process might be trapping or release of charges 
by impurity atoms or defects in the vicinity of the color cen- 
ters. Here, the number of trapped charges or varying distances 
of trap sites to the color center would affect the radiative 
quantum yield differently. The blinking time distribution of 
the Ni/Si centers shows remarkably long dark intervals of up 
to 500s. Such long time scales are known from investigations 
on photochromism of NV color centers l:26J and support the 
assumption of charge fluctuations. Anyway, the great length 
some dark periods renders a reliable usage of this color center 
as single photon source unfavorable. 
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Fig. 8 Saturation measurement of emitter (4). The background was 
recorded on the sample at a distance of 10|im from the color center. 



On the other hand, for emitter (4) we did not observe 
blinking but stable and bright emission. A measurement of 
the saturation count rate is displayed in Fig. |8] A fit for the 
count rate R of the form 



R^Ro. 



P + P., 



(8) 



to the data yields a saturation power Psat of 1.17mW and a 
saturation count rate R„ of 77.8kcount/s. This value is about 
50% lower than the 200kcounts/s reported for Ni/Si centers 
in nano-diamonds lfT2ll . which can be explained with the bet- 
ter optical coupling efficiency in the case of nano-crystalline 
diamonds. Regarding the focus size of our setup we can es- 
timate the saturation intensity as 365kW/cm^. For a rough 
estimate of the radiative excited state lifetime if of the emit- 
ter we assume a two-level system. Following Trebbia et al. 
Il27l T/ can then be calculated from Ti by using the relation 



(9) 



T/ = Ti- ( 1 + ) = 0.83ns -1.33= 1.11ns, 



which takes into account the excitation power of the corre- 
lation measurement (390|jW) and the saturation power (1.17 
mW). 

The lifetime of Ni/Si related centers in diamond nano- 
crystals on a silica substrate was reported to amount to 2 ns 
lfT2]| . The deviation of our result to this value can be explained 
with the radiative lifetime of a single emitter placed in a di- 
electric material. This lifetime decreases linearly with grow- 
ing refractive index n ll28l . In case of NV centers the radia- 
tive lifetime in bulk diamond (n = 2.4) was determined to 
be smaller by a factor of 2.2 than in nano-crystals on silica 
||29ll . We here find a similar factor of 1 .8 for the ratio of Ni/Si 
color center lifetime in bulk diamond and nano-crystals. For 
the precise determination of the time constants excited state 
lifetime measurements (via pulsed laser excitation) or excita- 
tion power dependent g''^' measurements are required which 
is work in progress. 



3.3 Conclusion 

We report the observation of single Ni/Si related color centers 
in single crystalline diamond produced via ion implantation, 
marking an important step towards their deterministic pro- 
duction. With different dose combinations for both species 
we found maximal doses of lO'^cm^^ for Ni and 10^" cm^^ 
for Si to yield a color center density that is as high as pos- 
sible and still keeps unwanted fluorescent background at a 
tolerable level. However, the creation efficiency is only about 
5 X 10^^ in relation to the Ni dose and 5 x 10^^ in relation to 
the Si dose and our results indicate that it does not depend on 
the implantation energy. The Ni/Si center emission lines are 
centered within a spectral interval from 767 to 775 nm with 
line widths between 1 .4 and 2.7 nm FWHM at room tempera- 
ture. Single photon emission could be demonstrated with val- 
ues for g'^' (0) that are clearly smaller than 0.5. The offset of 

gmLs (0) from zero can completely be explained with timing 
resolution of the HBT setup and the very short lifetimes of 
less than 1 .5 ns of the excited states of the emitters. 

The investigated emitters can be divided into two groups: 
The emitters of the first group show a double peak structure 
in the emission spectrum, photon bunching in the correla- 
tion measurement and a polarization visibility that is clearly 
smaller than 1 00 % in excitation. The color centers of the sec- 
ond group exhibit only one peak in the spectrum, they do not 
show any bunching and they have a polarization contrast of 
sa 98.5%. For the explanation of this different behavior we 
propose a three level system with excited and shelving state 
lying close to each other in energy. Strain in the diamond lat- 
tice or local charge distributions can possibly lead to close 
degeneracy, making thermalization possible and introducing 
an effective two-level system. 

For the Ni/Si centers we measured saturation count rates 
up to 78kcounts/s. Some of the centers showed strong three- 
level blinking which we tentatively explain by local charge 
fluctuations. 



4 Acknowledgements 

The ion implantations we were supported by Christian Hepp 
and Roland Albrecht (Universitat des Saarlandes, Saarbrii- 
cken, Germany) and performed with Prof. Dr Wolfgang Boi- 
se (Universitat Stuttgart, Germany) at the Universitat Gottin- 
gen, Germany and with PD Dr. Jan Meijer at the RUB ION 
(Ruhr-Universitat Bochum, Germany). The project was finan- 
cially supported by the Deutsche Forschungsgemeinschaft 
and the Bundesministerium fiir Bildung und Forschung (net- 
work EphQuaM, conti-act 01BL0903). 



References 

1 C. H. Bennett, G. Brassard, Proc. IEEE Int. Conf. on Computers, 
Systems & Signal Processing, 175 (Bangalore, India, 1984) 

2 A. Kuhn, M. Hennrich, G. Rempe, Phys. Rev. Lett. 89, 067901 
(2002) 



Single photon emitters based on Ni/Si related defects in single crystalline diamond 

3 C. Santori, M. Pelton, G. Solomon, Y. Dale, Y. Yamamoto, Phys. 
Rev. Lett. 86, 1502 (2001) 

4 A. M. Zaitsev, Optical Properties of Diamond (Springer, Berlin, 
2001) 

5 B. Lounis, M. Orrit, Rep. Prog. Phys. 68, 1 129 (2005) 

6 C. Kurtsiefer, S. Mayer, P. Zarda, H. Weinfurter, Phys. Rev. Lett. 
85, 290 (2000) 

7 F. Jelezko, T. Gaebel, I. Popa, M. Domhan, A. Gruber, J. Wrach- 
trup, Phys. Rev. Lett. 93, 130501 (2004) 

8 P. Neumann, N. Mizuochi, F. Rempp, P. Hemmer, H. Watanabe, 
S. Yamasaki, V. Jacques, T. Gaebel, F. Jelezko, J. Wrachtrup, 
Science 320, 1326(2008) 

9 V. A. Nadolinny, A. P. Yelisseyev, J. M. Baker, M. E. Newton, 
D. J. Twitchen, S. C. Lawson, O. P. Yuryeva, B. N. Feigelson, J. 
Phys.: Condens. Matter 11, 7357 (1999) 

10 T. Gaebel, I. Popa, A. Gruber, M. Domhan, F. Jelezko, J. Wrach- 
trup, New. J. Phys. 6, 98 (2004) 

HE. Wu, J. R. Rabeau, G. Roger, F. Treussart, H. Zeng, P. Grang- 
ier, S. Prawer, J.-F Roch, New. J. Phys. 9, 434 (2007) 

12 L Aharonovich, C. Zhou, A. Stacey, J. Orwa, S. Castelleta, 
D. Simpson, A. D. Greentree, F. Treussart, J.-F. Roch, S. Prawer, 
Phys. Rev. B 79, 235316 (2009) 

13 D. Steinmetz, E. Neu, C. Hepp, R. Albrecht, W. Boise, J. Meijer, 
D. Rogalla, C. Becher, Proc. SPIE 7727, 77270P (2010) 

14 C. Santori, P E. Barclay, K.-M. C. Fu, R. G. Beausoleil, Phys. 
Rev. B 79, 125313(2009) 

15 K.-M. C. Fu, C. Santori, P E. Barclay, R. G. Beausoleil, Appl. 
Phys. Lett. 96, 121907 (2010) 

16 C. Wang, C. Kurtsiefer, H. Weinfurter, B. Burchard, J. Phys. B 
39, 37 (2006) 

17 H. Sternschulte, K. Thonke, R. Sauer, P. C. Munzinger, P. Mich- 
ler, Phys. Rev. B 50, 14554 (1994) 

18 A. T. Collins, P M. Spear, J. Phys. C 16, 963 (1983) 

19 K. lakoubovskii, G. Davies, Phys. Rev. B 70, 245206 (2004) 

20 M. H. Nazare, A. J. Neves, G. Davies, Phys. Rev. B 43, 14196 
(1991) 

21 S. C. Kitson, P. Jonsson, J. G. Rarity, P. R. Tapster, Phys. Rev. A 
58, 620 (1998) 

22 I. Aharonovich, S. Castelletto, D. A. Simpson, A. Stacey, J. Mc- 
Callum, A. D. Greentree, S. Prawer, Nano Lett. 9, 3191 (2009) 

23 G. Davies, M. F Hamer, Proc. R. Soc. Lond. A 348, 285 (1976) 

24 A. T. Colhns, J. Phys.: Condens. Matter 1, 439 (1989) 

25 D. E. Gomez, J. van Embden, P. Mulvaney, M. J. Fernee, 
H. Rubinsztein-Dunlop, ACS Nano 3, 2281 (2009) 

26 K. lakoubovskii, G. J. Adriaenssens, M. Nesladek, J. Phys.: 
Condens. Matter 12, 189 (2000) 

27 J.-B. Trebbia, H. Ruf, P. Tamarat, B. Lounis, Opt. Express 17, 
23986 (2009) 

28 G. Nienhuis, C. T. J. Alkemade, Physica B+C 81, 181 (1976) 

29 A. Beveratos, R. Brouri, T Gacoin, J.-P. Poizat, P. Grangier, 
Phys. Rev. A 64, 061802 (2001) 



